Rocket dust storms and detached layers in the Martian atmosphere 
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Airborne dust is the main climatic agent in the 
Martian environment. Local dust storms play a 
key role in the dust cycle; yet their life cycle is 
poorly known. Here we use mesoscale modeling with 
radiatively- active transported dust to predict the evo- 
lution of a local dust storm monitored by OMEGA 
on board Mars Express. We show that the evolu- 
tion of this dust storm is governed by deep convec- 
tive motions. The supply of convective energy is 
provided by the absorption of incoming sunlight by 
dust particles, in lieu of latent heating in moist con- 
vection on Earth. We propose to use the termi- 
nology "rocket dust storm", or conio- cumulonimbus, 
to describe those storms in which rapid and effi- 
cient vertical transport takes place, injecting dust par- 
ticles at high altitudes in the Martian troposphere 
(30 to 50 km). Combined to horizontal transport by 
large-scale winds, rocket dust storms form detached 
layers of dust reminiscent of those observed with in- 
struments on board Mars Global Surveyor and Mars 
Reconnaissance Orbiter. Detached layers are stable 
over several days owing to nighttime sedimentation 
being unable to counteract daytime convective trans- 
port, and to the resupply of convective energy at sun- 
rise. The peak activity of rocket dust storms is ex- 
pected in low-latitude regions at clear season, which 
accounts for the high- altitude tropical dust maximum 
unveiled by Mars Climate Sounder. Our findings on 
dust-driven deep convection have strong implications 



for the Martian dust cycle, thermal structure, at- 
mospheric dynamics, cloud microphysics, chemistry, 
and robotic and human exploration. 



1 Introduction 

The Martian atmosphere has a permanent thin veil 
of suspended dust particles, the amount varying with 
location and season. The dust cycle is a crucial 
component of the Martian climatic system and the 



main factor of interannual variability (Smith 2004 



Montabone et al. 2005 ) . Airborne dust scatters and 



absorbs incoming sunlight and, to a lesser extent, out- 



going thermal radiation (Gierasch and Goody, 1972). 



This results in a significant warming of the Martian 
atmosphere, even in moderately dusty conditions. It 
is thus of key interest to characterize the spatial and 
temporal variability of dust in the Martian atmo- 
sphere and the control through lifting from the sur- 
face, transport by atmospheric winds and sedimenta- 
tion. 

The Martian tropospheric circulation is active at 
all spatial scales. This yields a variety of structures 
formed by dust lifted and transported by atmospheric 



winds (Cantor et al. 2001 Spiga and Lewis 2010 



Figure 1), which in turn change the Martian sur- 
face albedo ( Szwast et al. 2006 1 . Large-scale circula- 



tions, such as Hadley cells and baroclinic waves, in- 
duce dust transport over interhemispheric distances 
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(Wang et al. 2003). Mesoscale circulations, which 
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typical scales range from a few 100s kilometers to a 
few kilometers, give rise to dusty fronts along po- 



lar cap edges (Toigo et al. 2002), dust transport 



through slope winds (Rafkin et al. 2002) and local 



or regional dust storms (Malin et al.| 2008 Rafkin 



than GCMs. 

The interest for dynamical processes controlling 
the spatial distribution of dust in the Martian at- 
mosphere has been renewed by recent observations 
which completed pioneering observations on board 



2009 Maattanen et al. 2009). Turbulent eddies in Mariner and Viking (Anderson and Leovy 1978 



the Planetary Boundary Layer [PBL] cause gusts and 
vortices a few tens of meters wide, which can lift 
and transport dust, forming the frequently observed 



dust devils" (Thomas and Gierasch 1985 Balme 



and Greeley 20061. In addition to those phenom- 
ena encountered in the Martian atmosphere every 
year, a thick global dust loading occurs near peri- 



helion with irregular interannual variability (Zurek 



and Martin) [T9931 |Montabone et al.[ [20051 |Cantor 
2007). This phenomenon is often referred to as a 



"global dust storm" , despite the global dust cover not 
being related to an unique global storm but to the 
combined action of interhemispheric transport and 
regional dust storms growing unusually large. 

The Martian dust cycle is characterized by strong 
radiative-dynamical feedbacks. On the one hand, 
the amount of dust in the atmosphere is closely re- 
lated to atmospheric circulations through lifting and 
transport. On the other hand, airborne dust im- 
pacts thermal structure, hence atmospheric circu- 
lations at all scales. Thus far, the interplay be- 
tween dynamical and radiative phenomena control- 
ling the dust cycle have been mostly studied with 



Global Climate Models (GCMs) (Haberle et al. 1982 



Leovy 1985 ) . Dust-lifting parameterizations were in- 



cluded in GCMs to better explore those radiative- 



dynamical feedbacks (Murphy et al. 1995 Newman 



As pointed out by Rafkin ( 2009 ) , the GCM approach 



et al. 


2002 


Basu et al. 


2006 


Kahre et al. 


2006) 



suffers two main limitations. Firstly, GCM results 
are very sensitive to the choice of free parameters 
in the dust lifting and dust devil parameterizations. 
Secondly, lifting and transport phenomena within lo- 
cal, regional and global dust storms are associated 
with mesoscale circulations left unresolved by GCMs. 
Rafkin ( 2009 1 demonstrated through idealized three- 



dimensional simulations how mesoscale modeling per- 
mits to investigate the initiation and amplification 
of Martian dust storms through radiative-dynamical 
feedbacks in greater details and physical consistency 



Jaquin et al. 1986). Using Mars Express stellar oc- 
cupations, Montmessin et al. (2006) confirmed the 



occurrence of high-altitude detached hazes discovered 
in Mariner and Viking data, although the nature of 
aerosols (dust or ice) remained ambiguous. Mars Cli- 
mate Sounder [MCS] observations by McCleese et al. 
( 2010[ ) and |Heavens et al" (2011a) demonstrated that 
the vertical distribution of dust exhibits "detached 
layers" in apparent contradiction with the equilib- 
rium between large-scale mixing and sedimentation 



assumed thus far (Conrath 19751. Heavens et al 



(2011b I listed several possible scenarios (described in 
section 5.2 ) to account for those phenomena. In most 
cases, mesoscale and turbulent circulations appear 
to play a key role. Through Mars Global Surveyor 
[MGS] Thermal Emission Spectrometer [TES] mea- 
surements during the 2001 global dust storm, [Clancy| 



et al. (2010) also concluded that the observed dust 
vertical distribution implies "extraordinary vertical 
advection velocities" . 

In this context, further studies of dusty mesoscale 
circulations and their radiative-dynamical feedbacks 
are needed to broaden the knowledge of the Martian 
dust cycle and, more generally, of the Martian climate 
and meteorology. What is the evolution and fate of 
local dust storms on Mars? How do temperature and 
winds in the Martian atmosphere respond to such dis- 
turbances? What is their influence on the distribu- 
tion of dust particles? How could we characterize the 
mesoscale variability left unresolved by GCMs? Our 
original approach to address those key questions con- 
sists in modeling the evolution of a typical mesoscale 
storm observed in unprecedented detail by orbital in- 
frared spectrometry (Maattanen et al.l|2009|). A sum- 



mary of this observation is proposed in section[2] Our 
modeling strategy to simulate the storm is based on 
an upgraded version of the mesoscale model by[Spiga 
and Forget (2009), detailed in section [3] The pre- 



dicted evolution of the mesoscale dust storm and its 
impact on the dust distribution are analyzed in sec- 



2 



tion[4j We discuss the implications and perspectives 
of our findings in section [5] 



2 The OMEGA dust storm 

In a northern summer afternoon of Martian Year 27 
(L s = 135°, local time 1330), the OMEGA instru- 
ment on board Mars Express observed strong evi- 
dence of a local, optically very thick dust storm cen- 
tered at (25°E,2.5°S). Hereinafter we refer to this 
event as the "OMEGA storm" . The OMEGA storm 
took place in a region and at a season character- 



Visible dust optical depth 



ized by low dust storm activity (Cantor et al. 2001 



2006 Figure 17). Maattanen et al. (2009) used the 



OMEGA measurements to map optical thickness in- 
side the storm (Figure[T]) and to retrieve optical prop- 
erties of the Martian dust. An observation acquired 3 
sols earlier at the same location showed an essentially 
clear atmosphere, devoid of thick dust disturbances. 
The rapid onset and development of the OMEGA 
dust storm, peaking at visible optical thickness of 10 
in only three sols, are remarkable. Interestingly, using 
the same OMEGA observation three sols before the 
storm, Spiga et al. ( 2007 ) found strong perturbations 



in the surface pressure held retrieved from the 2 /xm 
CO2 line, suggesting that strong winds were already 
present and interacting with topography. 

Analysis of Figure [l] shows significant variability 
in optical depth within the OMEGA storm which 
exhibits a complex, convective-looking, mesoscale 
structure. Simultaneous Planetary Fourier Spec- 
trometer [PFS] measurements of temperature struc- 
ture also revealed how the Martian atmospheric state 
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Figure 1: Column dust optical depth at 1 /urn ob- 
served by OMEGA in Terra Meridiani. The zoomed 
area is located at coordinates (25°E,2.5°S). Values 
below the horizontal line in the colorbar cannot be 
considered as reliable. Adapted from Maattanen 



et al. (20091 with permission. 



is impacted by this dust disturbance (Maattanen Maattanen et 



detailed modeling study could help to reveal the evo- 
lution of the OMEGA storm. 

To complement OMEGA observations in 
((20091, we sought for MGS 



al. 



et al. 2009). Preliminary assessment of dust verti- Mars Orbiting Camera [MOC] images (Wang and 



cal distribution through single-column atmospheric 
modeling pointed towards low-altitude confinement 
of dust in the storm; yet these conclusions could 
not be fully reconciled with OMEGA nadir mea- 
surements. All these elements advocate for dedi- 
cated three-dimensional mesoscale modeling of the 
OMEGA storm. Moreover, further storm develop- 
ment and decay could not be followed by OMEGA 
in lack of Mars Express nadir data at the following 
ground track overlap three sols later. Thus, only a 



Ingersoll 2002 ) . Only global image mosaics acquired 



at local time 1400 could be used to attempt a 
tracking of the development and dissipation of the 
OMEGA storm. MOC images obtained for sols 
around the OMEGA storm are assembled in Fig- 
ure^ At L s — 132° the Martian atmosphere is clear 
in the region where the OMEGA storm occurred 
(marked with a red circle). At L s = 135°, close 
to the time of OMEGA measurements, the storm 
appears in the MOC image (unfortunately, MOC 
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images one sol before had a gap in observations 
exactly at the storm location). At L s — 136 — 137°, 
dust storm activity in the region continues during 
several sols and exhibits strong sol-to-sol variations 
and horizontal transport. This activity might be the 
result of either several individual short-lived (less 
than a sol) storms occurring in the same region, 
or one storm that was particularly variable both 
temporally and spatially. 



3 Modeling methodology 

3.1 Mesoscale model and general set- 
tings 

To study the evolution of the OMEGA dust storm, 
we employ the "Laboratoire de Meteorologie Dy- 
namique" (LMD) Martian Mesoscale Model ( MMM) 
(iSpiga and Forgetl [20091 ISpiga et al.l |2011|). This 



2009 


Spiga et al. 




2011) 


is built by coup. 


ing a 



mesoscale model 
dimensional, fully compressible, non-hydrostatic dy- 
namical core, capable to resolve fine-scale circula- 



tions on Earth (Skamarock and Klemp 20081, with 



the complete set of Martian physical parameteriza- 



tions used in the LMD-GCM (Forget et al. 1999 



Madeleine et al. 201U. The LMD-MMM has the 



capability to advect an arbitrary number of tracers. 
Initial and boundary conditions for the LMD-MMM 
are provided by LMD-GCM simulations which use 
similar physical parameterizations, thereby reducing 
inconsistencies. 

Details about the LMD-MMM and typical test sim- 



ulations can be found in Spiga and Forget (20091 



Physical parameterizations described in this refer- 
ence paper have been recently improved. The multi- 
size dust radiative transfer introduced by |Madeleine| 
et all (I2011I) is now employed in both the LMD-GCM 



and the LMD-MMM. Dust optical properties (extinc- 
tion function, single-scattering albedo and asymme- 
try factor) are computed using the most recent op- 
tical indices retrieved by Wolff et al. (2009). The 



advected tracers describing the distribution of dust 
aerosols are radiatively active, allowing for detailed 
analysis of the dust radiative-dynamical feedbacks. 
Additional details about this interactive dust scheme 








L = 136.73 
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Figure 2: Mars Daily Global Map retrieved with 
MOC instrument on board MGS. Images obtained 
with the "Mars Climate Center" website hosted by 
Ashima Research. 
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and the initialization of the OMEGA storm are pro- 



vided in sections 3.2 and 3.3 Multisize particle sedi- 
mentation is computed using Stokes' law and Cun- 
ningham correction factors (Rossow 1978). The 



shape of dust particles is accounted for by an addi- 
tional factor j3 in the second term of the Cunningham 
correction factor. We adopt j3 = 0.5 (disk-shaped 
dust particles) in the LMD-MMM as is done in the 
LMD-GCM (iMadeleine et all 120111). The evolution 



of the OMEGA storm predicted by the model is not 
sensitive to this parameter: similar results are ob- 
tained with j3 = 1 (spherical dust particles). 

The mesoscale domain employed for dust 
storm simulations is centered in Terra Meridi- 
ani at the approximate location of the OMEGA 
storm (25°E,2.5°S). It comprises 181 x 181 horizon- 
tal grid points with 7 km grid spacing (simulations 
with 10 km grid spacing were also performed and 
yield similar results) . Figure [3] shows the location 
and extent of the mesoscale domain projected on 
the Martian globe. Along the vertical dimension, 
101 levels are set with a spacing less than 1 km 
for altitudes below 40 km (Figure [3). LMD-MMM 
integrations are carried out with a timestep of 20 s, 
except for radiative transfer computations which 
are performed every 120 s (about l/31st of a 
Martian hour). This allows for radiative-dynamical 
feedbacks to be taken into account: changes in dust 
distribution rapidly impact the thermal structure. 
Mesoscale simulations are started at local time 1330 
and carried out in late northern summer (L s = 135°) 
when the OMEGA storm was monitored. Whenever 
different starting local time or L s are considered 



(section 4.2 and 5.1 ), initial and boundary conditions 
provided to the LMD-MMM by the LMD-GCM are 
recomputed. 

3.2 Interactive dust scheme and re- 
lated diagnostics 

The radiative impact at wavelength A of the total 
quantity of dust in the atmospheric column is given 
by the column optical depth t\ 



t\ -- I 'lr.\ 

I column 



(1) 
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Figure 3: LMD-MMM settings for dust storm sim- 
ulations. Top plot shows horizontal domain which 
comprises 181 x 181 horizontal grid points with 7 km 
grid spacing. Bottom plot shows vertical discretiza- 
tion which comprises 101 levels. Pressure at the top 
of the mesoscale domain is 1 Pa. x-axis denotes 
model vertical levels. j/-axis corresponds to (top left 
panel) mass-based coordinate, cf. Spiga and Forget 
( 2009[ ), (top right panel) altitude, (bottom left panel) 
pressure and (bottom right panel) distance between 
vertical levels. Values for altitude and pressure are 
computed with standard surface pressure 610 Pa and 
scale height 10 km; the actual model top in simula- 
tions is at ~ 55 km altitude. 
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The contribution dr\ to the column optical depth of 
a layer of dust of elementary thickness dP, where P is 
atmospheric pressure, writes 



dr\ = £ ^ cxt,A q dP with £ 



4pdust.9 



(2) 



where Pdust i s the dust particle density 
(about 2500 kg m -3 ), g the acceleration of gravity 
q the mass mixing ratio, r e g the effective radius, 
and Qext,A the dust extinction efficiency which is 
a function of r c g. For the sake of comparison with 
existing diagnostics in the literature, e.g. MCS 



retrievals (Heavens et al. 2011a I, the density-scaled 
optical depth 6 Z T\ may be employed instead of dr\ 

. -1 dr x gdr x Qext,\ ,„■. 

o z t\ = 1~ = ~T5~ = 9^ 1 3 ) 

p dz dP r e g 

where p is atmospheric density. In what follows, A 
subscripts are dropped for conciseness. Except oth- 
erwise indicated, dust optical depths are expressed 
at A = 0.67 /im. The conversion of dust optical 
depths and radiative properties from this wavelength 
to others, and vice-versa, is detailed in |Madeleine| 
et al. (2011). To compare quantities expressed 



at 0.67 /im with MCS observations in the A5 channel 
at 21.6 /im, we rely on the approximation 



7.3 



(4) 



as advised by Heavens et al. ( 2011a[) (this is consistent 
with Figure 1 in |Madeleine et al.[ |201l[ ). 



In LMD-MMM simulations, it is assumed that the 
sizes of the dust particles in any atmospheric layer 
follow a lognormal distribution of constant standard 



deviation ctq (Madeleine et al. 2011). Under this as 



sumption, the particle size distribution is fully de- 
scribed by the mass mixing ratio q and the dust num- 
ber density N (this method is often referred to as a 
"two-moment" scheme). The effective radius r c g is 
then proportional to the cube root of q/N 




with A 



(5) 



Thus the dust distribution in the LMD-MMM is 
merely predicted using two tracers [q, N], Those 



fields are modified at each model iteration by ad- 
vection (i.e. transport by winds), sources (lifting) 
and sinks (sedimentation). The radiative impact of 
the resulting dust distribution [q, N] is accounted 
for using equations [2] and [5] Through this in- 
teractive dust scheme, our mesoscale model offers 
a physically-consistent extrapolation of the storm 
structure and behavior from the instantaneous ob- 
servation by OMEGA. In what follows, the evolution 
of the simulated dust storm is diagnosed by the dust 
optical depth quantities t and S z t computed using 
equations [I] to [5] and predicted fields [q, N] in the 
LMD-MMM. 

3.3 Setting up the dust storm distur- 
bance 

LMD-MMM storm simulations are based on initial 
and boundary conditions obtained from LMD-GCM 
simulations. The interactive dust scheme described 
in section [3721 is used in both models. Column dust 
optical depths ttes measured by TES for a typical 
year devoid of global dust storms are employed to 
constrain [q, N] in the LMD-GCM, so that column 
optical depths Tb ac k in this model verify 



^back 5 



Q 



cxt 



column ' eft 



qdP = TTES 



(6) 



This dust distribution is thought to be representative 
of large-scale "background" dustiness. Initial back- 
ground fields in the LMD-MMM are denoted [q , N ] 
and T backi0 . 

To define an initial storm perturbation in the 
LMD-MMM, a local dust perturbation 5q is added to 
the initial background mass mixing ratio qo. The dust 
perturbation 5q is designed to mimic the column op- 
tical depth r s t orm of the mesoscale storm monitored 
by the OMEGA mapping spectrometer at 1 pm (Fig- 
ure [I]). It is defined such that 



Tstorm — £ 



Q 



cxt 



column r eS 



(q Q + Sq)dP (7) 



The initial dust perturbation 5q is assumed to ex- 
tend down to the surface where pressure is P sur f. We 
define Ptop as the pressure at the top of the dust 



6 



perturbation. Above the dust storm, for pressures P 
such that P < Ptop, we set Sq = 0. Within the dust 
storm, for pressures P such that Ptop < P < P SU rf> wc 
assume for simplicity that Sq is constant with pres- 
sure P, which yields 



Sq = (Tgtorm - Tback,o) 



surf Q. 



ext. 



top 



dP 



(8) 



The vertical distribution of dust within the 
OMEGA storm is not known. The rationale for 
choosing Sq constant with height within the storm is 
simplicity and the fact that afternoon PBL convec- 
tion yields well-mixed dust within a few hundreds of 
seconds (e.g. Spiga et al. 2010). Results are not criti- 



cally sensitive to this initial assumption. Conversely, 
the evolution of the OMEGA storm is sensitive to 
the pressure Ptop at the top of the dust perturba- 
tion, which is not constrained by OMEGA observa- 
tions either. The impact of the vertical extent of the 
storm disturbance is thus left to be explored in the 
modeling analysis (section [4.2.2 ). 

In the horizontal dimension, the storm area is a 
disk of center (25°E,2.5°S) and radius P s torm where 
the dust perturbation Sq is uniform. Outside the 
storm area, we set Sq = 0. To avoid sharp discon- 
tinuities at the storm boundaries, a transition zone 
of a few horizontal grid points is prescribed with hy- 
perbolic tangent functions between inside and outside 
the storm. 

Once the perturbed mixing ratio q' Q = qo + Sq is 
computed using equations [6] and |8j the number den- 
sity Nq follows from equation [5] We assume that the 
initial dust effective radius r e g is left unaffected by 
adding the storm perturbation upon background con- 
ditions. In the season and location considered for our 
mesoscale simulations, this yields r e g ~ 1.5 microns, 
which is also the value assumed by|Maattancn et al. 
(2009) for OMEGA retrievals. However, larger dust 



effective radius r e g are often observed in dust storms 
(IClancy et all 120101 lElteto and Toonl 120101). Hence 



we also carry out simulations where r e g is set to 3 fim 
within the storm perturbation (see section 4.2.2). In 
this case, dust optical properties, notably g cx t, are 
recomputed before equation [8] is used to obtain q' . 



4 Results 

4.1 Reference simulation 

The reference case study for the evolution of 
the OMEGA storm is defined as follows. The 
adopted dust optical depth in the storm perturbation 
is T s torm = 4.25 over an area of radius P s torm = 0.5° 
(~ 30 km). The value r s t orm ~ 4 does not ac- 
count for local maxima within the storm reaching 10, 
but reflects instead the averaged storm conditions 
monitored on board Mars Express by OMEGA (Fig- 
ure [T]) and PFS (see Maattanen et al. 



2009). We 



choose Ptop so that the initial dust storm perturba- 
tion is 10 km high, which corresponds to dust parti- 
cles being mixed within the convective PBL and the 
lowermost layers of the free troposphere. In the ref- 
erence simulation, dust lifting is not included: only 
the evolution of the storm perturbation witnessed by 
OMEGA is considered. The sensitivity of model pre- 
dictions to those various settings is discussed in sec- 
tion mu 

Figure [1] (left) shows the predicted column dust op- 
tical depth r every two hours from local time 1400 
to 0000. By the end of the afternoon, the dust storm 
appears to have drifted in the southwest direction, 
while its column optical depth r has been divided 
by a factor of two. What OMEGA saw as an op- 
tically thick dust storm in early afternoon becomes, 
according to LMD-MMM simulations, a dust distur- 
bance hardly discernible from background dust con- 
ditions in the middle of the night. In other words, our 
mesoscale simulations predict a rapid decay and short 
lifetime for the OMEGA dust storm. This tends to 
indicate that the spatial variability observed in MOC 
images in Figure [2] would correspond to several short- 
lived dust storms rather than a single, highly variable, 
dust storm event. 

The "nadir view" of the OMEGA dust storm in 
Figure [4] (left) does not reveal the most notable char- 
acteristics of its evolution. Figure [4] (right) shows 
the corresponding "limb view" through longitude- 
altitude sections of "MCS-like" density-scaled dust 
optical depth S z r mcs at 21.6 /im (cf. also movie 
in supplementary material). Within few hours, the 
dust storm disturbance evolves into a dusty plume 
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Figure 4: LMD-MMM reference simulation of the OMEGA storm, (left) Latitude-longitude maps of column 
dust optical depth t at 0.67 (im and (right) longitude-altitude sections of "MCS-like" density-scaled dust 
optical depth 6 z r mcs at 21.6 /im in 10~ 3 m 2 kg -1 . Sections are obtained at latitude 2.5°S. Season is late 
northern summer (L s = 135°). Diagnostics are shown every two hours from local time 1400 to 0000. 
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rising from lower to upper troposphere. The dusty 
plume reaches altitudes as high as 30 — 35 km around 
sunset. Then dust particles undergo significant hor- 
izontal transport and the storm loses its vertical 
extent to appear as a high-altitude elongated dust 
cloud. To summarize, mesoscale modeling predicts 
that within 6 hours the low-level dust disturbance 
monitored by OMEGA would turn into a fast-rising 
dust plume then a clear-cut high-altitude detached 
layer of dust. According to our model predictions, 
the implied horizontal and vertical advection of dust 
is spectacular: the vast majority of dust particles 
within the initial storm perturbation are transported 
well above the PBL to the higher troposphere and 
over horizontal distances reaching nearly a thousand 
kilometers (dust particles are even transported out 
of the mesoscale domain). This daytime dynamical 
transport accounts for the brief lifetime of the storm; 
sedimentation does not appear as a key factor. 

The plume-like appearance of the dust storm in 
Figure [4] denotes deep convective motions. This is 
confirmed by analyzing vertical winds predicted by 
the LMD-MMM. Figure [5] shows that updrafts within 
the dust storm often reach amplitudes > 3 m s _1 , 
peaking as high as 8 — 10 m s^ 1 , while compensating 
subsidences reach 1— 2 m s _1 . These values are about 
two orders of magnitude larger than typical vertical 
advection velocities obtained through GCM simula- 
tions (~ 10 _1 m s _1 ). Our predictions for vertical 
velocities in dust storms are in line with mesoscale 
simulations of |Rafkin ( |2009 ) (Figure 11) and Rafkin 



(2012) (Figure 3), as well as discussions in Clancy 



et al. (2010) and Heavens et al. (2011b) about pow- 



erful transport processes possibly accounting for the 
observed vertical distribution of dust on Mars. The 
predicted vertical wind field in Figure [5] also features, 
above the storm central updraft, alternating patterns 
reminiscent of upward-propagating gravity waves. 

Local dust storms such as the OMEGA storm have 
been sometimes referred to as "dust bombs" ; we pro- 
pose instead the terminology "rocket dust storms" to 
better emphasize how rapid and powerful the verti- 
cal transport can be. Alternatively those phenomena 
can be named conio-cumulonimbus, using konios the 
Greek word for dust (Heavens et al. 2011b|). This 




East Longitude <deg) 



emphasizes the analogy with pyro-cumulonimbus on 



Figure 5: LMD-MMM reference simulation of the 
OMEGA storm. Longitude-altitude sections ob- 
tained at local time 1500 and latitude 2.5°S. From 
top to bottom: (a) MCS-like density-scaled optical 
depth at 21.6 fim in 10 -3 m 2 kg -1 ; (b) Vertical wind 
in m s _1 (maximum is about 10 m s _1 ); (c) Short- 
wave heating rate in K per Martian hour (maximum 
is about 24 K/h). 
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Earth, associated with wildfires and thought to be 
responsible for detached layers of aerosols in the ter- 



restrial stratosphere (Fromm et al. 2010) 



Absorption of incoming sunlight by transported 
dust is the driving mechanism for the deep con- 
vective motions in rocket dust storms, playing the 
role devoted to latent heat in terrestrial moist con- 
vective storms. Radiative heating rates %sw hi 
shortwave (visible) wavelengths are displayed in Fig- 
ure [5] Shortwave absorption by dust particles trans- 
ported within the storm yields an extreme warming 
of 15 — 20 K per Martian hour in atmospheric lay- 
ers where density-scaled optical depth is large. Be- 
low those layers, and at the surface, extinction yields 
moderate shortwave radiative warming, which does 
not appear to lead to significant dynamical effects. 
Absorption by dust particles in infrared wavelengths 
also warms the atmosphere though one order of mag- 
nitude less than shortwave absorption. Radiative 
warming within the dust storm disturbance causes 
the atmosphere to become positively buoyant, which 
gives rise to strong upward convective motions. 

Figure [6] shows through a latitude-altitude section 
how rocket dust storms impact thermal structure. 
Compared to background temperature conditions, at- 
mosphere is warmer above/within the maximum in 
density-scaled optical depth and cooler below it. This 
appears in qualitative agreement with PFS tempera- 
ture measurements over the OMEGA storm (Figure 
8 in Maattanen et al. 20091 which features a posi- 



tive temperature anomaly in the troposphere and a 
negative one near the surface. This probably means 
that the OMEGA storm has already started its as- 
cent at the initial local time considered in reference 



simulation (this scenario is explored in section 4.2.2 



and does not change the results presented here). 

The impact of rocket dust storm on temperature 
structure is best accounted for by analyzing poten- 
tial temperature in Figure [6j From the Lagrangian 
point of view (following an ascending dusty parcel) , 
potential temperature 9 is only modified through di- 
abatic processes J 



D6 
~Dt 



= J 



with 



= TV' 1 and J = 



H 



(9) 
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Figure 6: LMD-MMM reference simulation of the 
OMEGA storm. Latitude-altitude sections obtained 
at local time 15 : 00 and longitude 25°E. (top 
plot, shaded) MCS-like density-scaled optical depth 
in 10~ 3 m 2 kg" 1 , (top plot, contours) Temperature 
in K. (bottom plot, shaded) Vertical wind in m s _1 . 
(bottom plot, contours) Potential temperature in K. 
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where c p is specific heat capacity, T temperature, 
t time, 1-L diabatic heating rate, V = (P/ Pq) R ^ Cp 
Exner function, R gas constant, Pq reference pres- 
sure. This formulates how diabatic warming turns 
the dust disturbance into a warm bubble of poten- 
tial temperature which subsequently rises. Here, 
H — Hsw and variations of "Hsw mostly follow vari- 
ations of density-scaled optical depth (see Figure [5]) . 
Now, from the eulerian point of view, potential tem- 
perature at a given geometrical position in Figure [6] 
is modified both through diabatic processes J and 
vertical advection of heat A 

dO dO 

— = J + A with A = -w— (10) 

at oz 

where w is vertical wind (horizontal transport of 
heat is neglected) and z altitude. At the top of the 
dust disturbance in Figure |6j potential temperature 
is larger than in the environment: J is positive owing 
to radiative warming by dust particles while A is close 
to zero because convective motions are not yet estab- 
lished. Conversely, at the bottom of the dust dis- 
turbance, potential temperature is lower than in the 
environment: A is negative, since established vertical 
winds induce advection of air with lower from be- 
low, and overwhelms J as dust particles (drivers for 
diabatic heating J7) are being transported upward. 
Hence the temperature and potential temperature 
disturbances in Figure [6] are governed by competing 
radiative warming and adiabatic cooling through as- 
cent. Only fine-resolution mesoscale modeling is able 
to resolve this interplay between radiation and dy- 
namics. We also note that emitted gravity waves im- 
pact temperature fields above the rocket dust storm. 

In other words, rocket dust storms are radiatively- 
controlled dry convective storms. Similarly to ter- 



restrial moist convection (see e.g. Holton 20041, we 



shall define the Convective Available Potential En- 
ergy [CAPE] C associated with Martian dust storms 



C = 



AT , 

9 ~ — dz 



T 

storm - 1 env 



(11) 



where T env is the environmental temperature (out- 
side the storm) and AT is the temperature contrast 
between inside and outside the storm disturbance (i.e. 



the temperature of rising air parcels minus the en- 
vironmental temperature). Assuming ideal energy 
transfer and negligible entrainmcnt, the maximum 
vertical velocity ui max reached by convective updrafts 
writes 

w max = V2C (12) 

A 5 km-deep dust plume undergoing, as in Figure [5] a 
radiative warming of AT ~ 7 K during half an hour 
while its environment is at 250 K would rise verti- 
cally with typical velocity u> m ax ~ 30 m s _1 (cor- 
responding CAPE is C ~ 500 J kg" 1 ). Clearly this 
theoretical calculation is an overestimate, in which 
significant entrainment and rapid plume rising are 
neglected, but the order of magnitude obtained from 
LMD-MMM simulations in Figure [5] is correctly ac- 
counted for. The meaning of this simple calculation 
is qualitative rather than quantitative: it illustrates 
that dust radiative effects have a strong potential to 
trigger deep convection on Mars. 

Apart from large temperature contrast between 
inside and outside storm, rocket dust storms could 
reach high altitudes for two main reasons. Firstly, the 
upward transport of dust particles (and subsequent 
rising of the area where radiatively-induced warming 
is maximum) acts as a positive feedback for storm 
convective motions. Mesoscale simulations show that 
induced shortwave heating rates does not vary much 
when the OMEGA storm rises. Since environmental 
temperature decreases with altitude, this results in 
larger C hence faster updrafts. Secondly, processes 
which could inhibit convection are not efficient. At 
pressures as low as 10 Pa, Stokes-Cunnigham sedi- 
mentation velocities for dust micron-sized particles 
are still one to two orders of magnitude lower than 
convective velocity. Moreover, on Mars, no stable 
layer like the stratosphere on Earth is encountered 
by rising dust plumes whose summits are, as a con- 
sequence, less likely to resemble the "anvil" of terres- 
trial cumulonimbus. What really causes the rocket 
dust storm to stop rising is the abrupt drop in in- 
coming sunlight at sunset. 

Despite strong vertical acceleration, the OMEGA 
storm is still impacted by horizontal winds as is 
shown in Figure [4] (and movie in supplementary ma- 
terial). This is in line with MOC images in Figure [2] 
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Figure 7: LMD-MMM reference simulation of the OMEGA storm. Maps of horizontal wind speed in m s _1 
with wind vectors superimposed and density-scaled optical depth contoured. Top plots show local time 1400 
and bottom plots local time 1500. Altitude above MOLA zero datum is 3 km (top left), 14 km (top right), 
10 km (bottom left), 20 km (bottom right). 
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Tropospheric jet-streams act to gradually drift dust 
particles towards southwest, before the storm really 
loses its plume-like appearance after sunset and the 
sudden drop in convective energy supply. Mesoscale 
temperature gradients induced by the dust storm it- 
self contribute to horizontal transport of dust parti- 
cles. Hydrostatic equilibrium implies that those tem- 
perature gradients drive outward (inward) pressure 
forces at the summit (bottom) of the storm, there- 
fore divergent (convergent) horizontal flow in the ab- 
sence of significant Coriolis force. LMD-MMM pre- 
dictions in Figure [7] show that horizontal inflow and 
outflow induced by rocket dust storms are significant 
and yield local modulations of large-scale winds by 
about a factor of 2. These thermally-induced hori- 
zontal motions cause a confinement of dust particles 
in the lower part of the storm and a dust storm widen- 
ing near the summit. 

Figure [8] displays the evolution of the dust cloud 
in the night and in the morning the following sol; 
maps of dust optical depth are omitted here because, 
from the nadir point of view, the high-altitude dust 
cloud is no longer discernible from background con- 
ditions. Figure [8] shows that dust particles initially 
composing the OMEGA storm in the afternoon form 
a clear-cut detached layer at nighttime. This de- 
tached layer of dust is predicted to occur at alti- 
tudes 25 — 30 km and to extend over large horizontal 
distances (more than 5° degrees longitude). It ex- 
hibits structure and altitude reminiscent of detached 
layers observed by MCS at pressure 10 — 100 Pa 
(Heavens et al. 2011a|b I. Moreover, typical maxi- 
mum values for density-scaled optical depth <5 2 r mcs 
in detached layers is predicted by LMD-MMM sim- 



ulations to reach ~ 3 to 6 x 10 m kg 



fair 



agreement with MCS observations of enriched lay- 
ers of dust which range from 1 to 4 x 10 -3 m 2 kg -1 
[cf. 



Figure 
in 



4 in 



ure 7 
servations 



Heavens et al 



(|2011aj) and Fig- 
Heavens et al.| (|2011bj)]. Both MCS ob- 
and mesoscale simulations indicate that 
density-scaled optical depth within high-altitude en- 
riched layers of dust is about 2 to 10 times the value in 
lower troposphere. We thus show that dust storms in 
the lower troposphere, through their "rocket storm" 
behavior in the afternoon, have the ability to lead 
to detached layers of dust in the upper troposphere 




Figure 8: Same as Figure [4] except that section lat- 
itude is 4°S and local times 02 : 00 to 12 : 00 every 
two hours (from top to bottom). 
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within a few hours. 

Furthermore, LMD-MMM simulations give in- 
sights into the lifetime of detached layers of dust. 
Figure [8] shows that the predicted high-altitude en- 
richment in dust particles in the morning is only 
slightly lower than in the previous evening. The 
most noticeable change of this dust cloud in six 
hours is instead a 5-km descent through sedimenta- 
tion. This corresponds to downward vertical velocity 
of 0.2 m s _1 , consistent with the above-mentioned 
Stokes-Cunnigham estimates. Enhanced radiative 
cooling yields a reversed thermal circulation in the 
dust cloud, which moderately reinforces downward 
motions in nighttime. This scenario for nighttime 
evolution could explain why dust detached layers 
are often located at higher altitude in daytime than 



in nighttime in MCS observations (Heavens et al. 



2011a, though the authors recognize observational 
challenges might cast shadow on this conclusion). 

There is a strong asymmetry between mesoscale 
processes leading to vertical transport of dust parti- 
cles. Through convection, it only takes few hours for 
dust particles to rise from the lowermost troposphere 
to altitudes as high as 35 km; through sedimentation, 
it would take at least several sols for them to come 
back to their initial altitude. A detached layer of dust 
which persists in the Martian atmosphere the whole 
night is actually likely to stay until the following 
evening (provided no sudden increase in large-scale 
winds is encountered). Figure [8] shows indeed that, 
after sunrise, dust detached layers exhibit convective, 
plume-like, structures caused by the resupply of con- 
vective energy through absorption of visible sunlight 
by dust particles. This causes the dust detached layer 
in our case study to rise to altitudes 30 — 40 km. 
Through this mechanism, high-altitude enriched lay- 
ers of dust are expected to survive for at least several 
sols, if not tens of sols, before changes in large-scale 
winds would cause density-scaled optical depth to 
drop as horizontal winds spread dust particles over 
large distances. There resides the potentially long 
lifetime of detached layers of dust in the Martian at- 
mosphere. 



4.2 Sensitivity simulations 

Figure [9] summarizes the evolution of the OMEGA 
storm predicted by the LMD-MMM with reference 
settings: rocket dust storm in the afternoon, for- 
mation of a detached layer in the evening, descent 
through sedimention in the night, ascent in the fol- 
lowing morning which ensures detached layers are 
able to survive during several days. The robustness 
of this dynamical scenario is now confronted to addi- 
tional simulations with modified settings. 



4.2.1 Radiative forcing of storm disturbance 

Absorption of incoming sunlight by dust particles is 
the supply of convective energy in rocket dust storms. 
A key sensitivity test is to assess the OMEGA storm 
evolution without the radiative effects of transported 
dust. To that purpose, we keep the column dust op- 
tical depth equal to TES column optical depth (i.e. 
background dust optical depth) in LMD-MMM radia- 
tive transfer calculations. Initial storm disturbance 
and LMD-MMM integrations are similar to the ref- 
erence simulation. Dust particles in the storm per- 
turbation are transported by atmospheric winds and 
undergo sedimentation, but play no radiative role. 

From the nadir point of view, the storm lifetime is 
similar whether or not transported dust is radiatively 
active (in the two cases, column optical depth looks 
as in Figure [4]) . Storm lifetime is slightly longer in 
the radiatively-active case, owing to induced thermal 
circulations which imply convergence at the bottom 
of the storm and limit the horizontal spreading of 
dust particles by large-scale winds. Conversely, the 
limb point of view reveals major differences between 
the radiatively active and inactive cases (compare 
Figure 10 with Figure HI) . Deep convective motions 



(rocket dust storm) and subsequent high-altitude en- 
riched layers of dust are not predicted when trans- 
ported dust is assumed to be radiatively inactive. In 
the absence of strong vertical winds, horizontal wind 
shear dominates the transport of dust particles. The 
initial dust disturbance is almost entirely dissipated 
in the end of the night under the influence of hori- 
zontal transport and sedimentation. Radiative effects 
of transported dust appear essential to initiate and 
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Figure 9: LMD-MMM reference simulation of the OMEGA storm. This figure provides a summary of the 
storm evolution: rocket dust storm in the afternoon (top left panel), formation of a detached layer in the 
evening (top right panel), descent through sedimention in the night (bottom left panel), resuming ascent 
in the following morning which ensures detached layers are able to survive during several days (bottom 
right panel). Longitude-altitude sections of "MCS-like" density-scaled dust optical depth are shown as in 
Figures [4] and [8j except units are 10~ 4 m 2 kg -1 and results are averaged over latitudes (0 — 5°S) so that the 
section latitude has not to be adapted to the southward drift of the OMEGA storm. The displayed local 
times are also distinct from Figure [4] and indicated on each panel. 
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maintain high-altitude detached layers of dust in the 
Martian atmosphere. 



4.2.2 Properties of initial storm disturbance 




Here we assess how sensitive the predicted evolution 
of the OMEGA storm is to assumptions made on the 
initial dust disturbance. The sensitivity simulations 
carried out to that purpose are listed in table [TJ The 
reference case is case R. In cases n and nD, a more 
confined dust disturbance is considered to simulate 
the observed local maximum in dust optical depth in 
Figure [l] Case D is similar to case R except that 
the storm optical depth is larger. Case W explores 
the possibility that the OMEGA storm might actu- 
ally be larger than the OMEGA footprint. Case LP 
is an alternative simulation designed to account for 
larger dust effective radius r e g often observed in dust 
storms (see section 3.3). Case IT assumes an ini- 
tial 5 km-high dust storm perturbation instead of 
the reference 10 km-high perturbation. Indeed, only 
transport and mixing below PBL top (~ 5 km in the 
considered location and season) can be undoubtedly 
assumed: assuming that the initial dust disturbance 
extends above the PBL top is implicitely assuming 



Figure 
within 



that deep convective mixing occurs (see section 4.1). 
Case E considers that the storm disturbance is not 
confined near the surface at the local time of OMEGA 
measurements (1330): as is suggested by both the 
analysis in section |4.1| and PFS measurements, the 
convective ascent of the storm has probably already 
started at that local time. Given the storm convec- 
tive behavior, this scenario is equivalent to setting a 
disturbance similar to the reference case at an earlier 
local time (we chose 1130). 
»» L °™«" d '»«> The results of sensitivity simulations are shown in 

Figure |11| for comparison with the reference case in 
Figure |9j (results for case D are close to case nD hence 

10: Same as Figure HI except dust particles omitted )- For a11 cases in Table the sto ™ evolu " 
storm are not radiatively active. tion is qualitatively similar to the reference case: af- 

ternoon convection (rocket dust storm), evening de- 
tached layer, nighttime sedimentation, morning re- 
ascent which acts to maintain the detached layer. 
The reference simulation appears as typical of what 
could be obtained for a range of reasonable assump- 
tions for the initial dust disturbance. A few elements 
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Case 


Comments 


Tstorm 


^?storm 


r e & 


Bottom-top altitudes 


R 


Reference 


4.25 


0.5° (~ 


30 km) 


1.5 fj,m 


- 10 km 


n 


Narrow 


4.25 


0.1° (~ 


; 6 km) 


1.5 fj,m 


- 10 km 


nD 


Narrow dustier 


10 


0.1° (~ 


'. 6 km) 


1.5 //m 


- 10 km 


D 


Dustier 


10 


0.5° (~ 


30 km) 


1.5 /jm 


- 10 km 


W 


Wide 


4.25 


1.0° (~ 


60 km) 


1.5 /im 


- 10 km 


LP 


Larger particles 


4.25 


0.5° (~ 


30 km) 


3.0 /im 


- 10 km 


IT 


Lower top 


4.25 


0.5° (~ 


30 km) 


1.5 /jm 


- 5 km 


E 


Earlier storm 


4.25 


0.5° (~ 


30 km) 


1.5 fim 


- 10 km [1130LT] 



Table 1: Parameters defining the initial dust disturbance and explored in the sensitivity study detailed in 
section [472] and Figure [TT] 






Case LP Case IT Case E 




Figure 11: Figures analogous to Figure [9] for cases detailed in table [T] and described in section [472] Note 
that the displayed local times for case E are earlier than other cases (hence morning re-ascent is not shown 
contrary to other displayed cases). In cases n and nD (case W) lower (larger) values for dust-scaled optical 
depth compared to Figure [9] are caused by meridional averaging since storm radius is lower (larger) . In 
cases LP and E, density-scaled optical depth in the morning dust cloud is less than in the reference simulation 
because the rocket dust storm reached higher altitudes characterized by enhanced horizontal winds which 
advect more dust particles out of the mesoscale domain. 
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in Figure 11 are worth being detailed. 



• The size of initial perturbation does not impact 
the altitude reached by OMEGA dust storm by 
the end of the afternoon (cases nD and W) . The 
main difference with the reference case R is the 
quantity of dust particles injected at high alti- 
tude (also worthy of notice in case W, 2 — 3 dis- 
tinct rising plumes develop instead of one) . 

• A dustier storm injects particles at higher alti- 
tudes (cases D and nD). This is expected since 
radiative warming hence supply of convective en- 
ergy are larger. The altitudes reached by dust 
storm particles in those cases (45 — 50 km) are 
probably overestimated though. Figure [I] indi- 
cates that the r s torm = 10 area is embedded into 
an area where optical depth is larger than as- 
sumed in cases nD and D where temperature 
contrast AT with the environment, and resulting 
CAPE, are therefore overestimated. 

• The altitude of detached layers of dust in the 
night does not vary much with the initial dusti- 
ness of storm perturbation (compare e.g. cases n 
with nD). Dust particles reach higher altitudes 
in the end of the afternoon, but lower environ- 
mental densities cause larger sedimentation ve- 
locities in the night. 

• The high-altitude enrichment of dust in the 
end of the afternoon is enhanced and occurs 
at higher altitudes if dust effective radius is 
assumed to be larger (case LP). Computation 
of dust radiative properties as in | Madeleine" 



et al. (2011 1 shows that larger r e fj yields slightly 



lower shortwave Qext/'cff ratio and larger long- 
wave Qext/feS ratio. Heating rates follow simi- 
lar trends (equation^. LMD-MMM simulations 
show indeed that increasing r e g from 1.5 /im 
to 3 fiui raises longwave heating rates from ~ 2 
to ~ 8 K per Martian hour, while variations 
in shortwave heating rates are marginal (less 
than 10%). This yields enhanced CAPE C within 
the rocket dust storm, hence enhanced vertical 
transport of dust particles. This is however mit- 
igated in nighttime by enhanced sedimentation 



rates associated with larger r e g, as well as a more 
efficient thermally-induced circulation caused by 
enhanced longwave cooling. As a result, the al- 
titude of dust detached layer in the end of the 
night is roughly similar to case R. 

• The altitude reached by the rocket dust storm 
is 10 km lower for an initial perturbation more 
confined near the surface (case IT) than the ref- 
erence case R. This apparently contradicts the 
larger supply of convective energy which stems 
from the larger density-scaled optical depths 
associated with near-surface confinement (pre- 
dicted heating rates reach 30 K per Martian 
hour). Nevertheless, enhanced thermal circula- 
tions are associated with the warmer dusty dis- 
turbance, which causes reinforced divergence at 
storm top hence horizontal outflow of dust par- 
ticles. As a consequence, the detached layer of 
dust in case IT forms at 20 — 25 km, at lower 
altitude than case R, but over larger latitudinal 
extent. 

• When the storm simulation is started at local 
time 1130 (case E), the dust disturbance has 
risen to pressure ~ 100 Pa at local time 1330, 
in agreement with PFS measurements (cf. the 
approximate altitude of transition between neg- 
ative and positive temperature anomalies in Fig- 
ure 8 of Maattanen et al. (2009)). Since the 



ascent of the rocket dust storm begins earlier 
than in the reference case, the maximum altitude 
reached is larger. Apart from this difference, the 
evolution of the dust detached layer in case E is 
similar to the reference case R. 

4.2.3 Dust lifting 

Our study focusses on the evolution of an observed, 
already established, dust storm perturbation. It is 
out of the scope of this paper to assess how the 
OMEGA storm initially appeared. Triggering a phe- 
nomenon such as the OMEGA dust storm probably 
requires a peculiar combination of large-scale and 
mesoscale processes. It is challenging to reproduce 
this with Martian mesoscale models which have not 
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yet reached the accuracy of their terrestrial counter- 
parts. Furthermore, mechanisms governing the injec- 
tion of Martian dust particles beyond the lowest cou- 
ple of meters in PBL is still an open problem (cf. re- 



and Bergametti ( 1995 1 



views by Greeley and Iversen 1987 Kok et al. 2012 1 



Notwithstanding this, a less idealized case than 
those considered in sections 14.11 and 14.21 can be de- 
signed to assess both the birth and evolution of the 
OMEGA storm. The LMD-MMM is run in the same 
conditions as in the reference case, except no initial 
dust perturbation is assumed: instead dust lifting 
through near-surface winds is activated over the area 
covered initially by the OMEGA storm perturbation 
in the reference case. This allows us to simulate a 
more realistic evolution of the OMEGA storm (no- 
tably, possible feedbacks in which storm circulations 



impact dust lifting as in Rafkin 2009) without ad- 



dressing the complexity of storm initiation which re- 
quires dedicated studies. 

A simple lifting scheme is adopted. The lifting of 
dust particles by near-surface winds is mainly a func- 
tion of wind stress a = pul exerted on the surface. 
Friction velocity it* can be approximately deduced 
from horizontal wind speed u(z\) predicted a few me- 
ters above the surface 



ku(zi) 



(13) 



where k = 0.4 is von Karman's constant and z$ the 
roughness length (set to 1 cm in the LMD-MMM). To 
first order, the threshold stress at = pv% t (with u*t 
the threshold friction velocity) required to initiate 
particle motion depends on particle size. Models 
and wind tunnel experiments show that near-surface 
winds tend to catch large dust particles (10-100 fxm) 
rather than micron-sized dust particles suspended 



in the Martian atmosphere (White 1979 Newman 



et al. 2002|. Thus saltation and sand-blasting are 



possibly playing a key role in dust lifting on Mars 



( Greeley 2002 ) . Hence the simplest estimate of at for 
micron-sized particles is the value for large particles, 
which approximately ranges between 10 mN m 
and 40 mN m -2 . Once this threshold is reached, 
the vertical flux of dust particles V^ft (kg m -2 s _1 ) 



P 2 

Vuft = 2.61 a - (tt* - it*t) (u* + u*t) (14) 
9 



where a 
range [10 



is an 

-Mo- 1 ! 



efficiency coefficient in the 
m . If a > at, dust mass 
mixing ratio q is modified in lowermost model layers 
given mixing coefficients and the vertical flux Vuft of 
dust particles lifted from the surface. The increase 
in number density N is obtained assuming lifted 
dust particles follow a lognormal distribution with 



effective radius 3 /im (see sections 3.3 and 4.2.21. 

The values of lifting threshold at and efficiency a 
appropriate for Mars are unknown and certainly not 
uniform over the planet. Yet the detection of the 
OMEGA storm hints at favorable lifting conditions 
(owing to either availability of dust at the surface 
or exceptional meteorological conditions) . Hence low 
threshold is assumed: at — 5 mN m~ 2 . It is cho- 
sen so that u*t ~ 0.5 m s _1 , which ensures predicted 
peaks in near-surface winds imply lifting of dust par- 
ticles from the surface (this optimistic value might 
not b e so if hysteresis effects are at play, see Kok[ 



2010). We set a = 2 x 10 m inside storm area 



(a typical value according to literature, Kahre et al. 
20061 IRafkinl 120091) and a = outside. First-order 



calculations with equation 14 and chosen [a, at] yield 



vertical flux Vjift compatible with dust mass mixing 



ratio q within storm calculated in section 3.3 



The LMD-MMM simulation is started in early af- 
ternoon as in the reference case. Dust lifting starts 
during the night under the influence of slope winds 
associated with cratered terrains. Contrary to day- 
time conditions, dust particles are confined near the 
surface in the night: ultra-stable conditions lead to 
shear-driven PBL mixing with limited vertical extent. 
This is analogous to the adverse influence of temper- 
ature inversion layers on the vertical transport of pol- 



is estimated following White ( 1979 ) and Marticorena 



lutants on Earth. Figure 12 shows that this results in 
a significant enrichment of dust in the lowermost PBL 
layers, the optical depth thereof reaching r ~ 3. Hor- 
izontal transport of dust is northward in early morn- 
ing but turns southwestward towards the afternoon. 
During the morning, as the sun rises in the Mar- 
tian sky, dust particles are being transported upward 
increasingly quickly as the near-surface dust distur- 



19 



bance turns into a rocket dust storm (Figure 12). In- 
creasing incoming solar flux enhances the supply of 
convective energy through absorption by dust parti- 
cles initially at low altitudes. Later in the afternoon, 
the resulting extinction of incoming sunlight passing 
through the top dust-rich layers does not prevent suf- 
ficient radiative warming to occur below. The day- 
time development of the rocket dust storm is also 
facilitated by the growth in PBL mixing depth up to 
about 5 — 6 km above the surface. 

The convective behaviour of dust perturbations is 
associated with low-level convergence and enhance- 
ment of near-surface winds which in turn induce lift- 
ing of dust particles from the surface. This is reminis- 
cent of the positive radiative-dynamical feedbacks de- 



scribed in Rafkin (2009), possibly giving rise in mid- 



latitudes to balanced circulations analogous to ter- 
restrial hurricanes. In our mesoscale simulations, the 
feedback mechanism is however limited because the 
fast ascent of rocket dust storms yields low duration 
of storm-induced near-surface circulations capable to 
lift dust. A balanced low-level dust storm cannot es- 
tablish in equatorial and tropical regions, as is also 
noted by Rafkin ( 2009 ) . This is confirmed in our sim- 



ulations by the absence of low-pressure cores in rocket 
dust storms. The short lifetime of the OMEGA storm 
is not caused by inefficient lifting (given our chosen 
threshold at) but instead by the fast, convective, as- 
cent of the dust disturbance. 

Dust particles are being transported upward to al- 
titudes 25 — 35 km within 2 — 4 hours to eventu- 
ally form a dust detached layer. Figure [12] can be 
compared to Figure g] The evolut ion of the storm 
disturbance is qualitatively similar to the reference 
case, except that the storm ascent seems to have al- 
ready begun at local time 1330, which further jus- 
tifies the sensitivity simulations carried out in sec- 
Predictions in Fig- 



tion 
ure 



4.2.2 (case E in Figure 11) 



12 after a full simulated Martian sol (and simi- 



lar predictions obtained in third simulated sol) also 
indicate that LMD-MMM results are not too sensi- 
tive to spin-up (i.e. the few hours at the beginning 
of a mesoscale simulation during which predictions 
adjust to GCM boundary conditions). All these el- 
ements tend to confirm the robustness of the storm 
evolution derived from the reference case, despite the 



idealized character of its initial state. 

5 Discussion 

In the particular case of the OMEGA storm, 
LMD-MMM simulations demonstrate that local dust 
storms have the ability to form detached layer of 
dust through radiatively-induced convective motions, 
a behavior we proposed to name "rocket dust storm" . 
Physical mechanisms discussed for this particular 
case are general enough for similar rocket dust storms 
to appear in other regions and seasons. In this sec- 
tion, we discuss the expected spatial and seasonal 
variability of rocket dust storms, their impact on the 
vertical distribution of dust (notably the formation 
of detached layers of dust) and other consequences 
on the Martian dust cycle, thermal structure, atmo- 
spheric dynamics, cloud microphysics, chemistry, and 
robotic and human exploration. 

5.1 Variability of rocket dust storms 
and subsequent detached layers of 
dust 

Where and when are rocket dust storms expected 
to appear? For reasons detailed in section 4.2.3 we 



discuss other factors than the necessary condition of 
dust being available for lifting and lifting threshold 
being reached. For similar reasons, we consider here 
a given storm disturbance and do not discuss vari- 
ability caused by the quantity of lifted dust particles 
initially composing the storm disturbance. 

In rocket dust storms, convective energy is pro- 
vided through absorption of incoming solar radiation 
by suspended dust particles. This makes solar energy 
input a crucial factor governing rocket dust storms: 
the further the dust disturbance from subsolar lati- 
tudes, the weaker the convection, the lower the alti- 
tude reached by transported dust particles. 

A second key factor to control the strength of 
rocket dust storms is background dustiness. Con- 
vective motions would be particularly developed for 
large CAPE C, which happens for large temperature 
contrasts AT between the storm and its environment. 
For a given input in solar energy, the less dusty the 
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Figure 12: LMD-MMM storm simulation with lifting and no initial dust perturbation. Same as Figure [4] 
except that local times range from 0800 to 1800 and longitude-altitude sections are obtained at latitude 1.5°S. 
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background atmosphere, the lower the environmental 
temperature, the larger the vertical winds and alti- 
tudes reached by rocket dust storm. The OMEGA 
storm occurred at a season (late northern summer) 
when the Martian atmosphere is clear which, in ad- 
dition to large incoming solar flux, explains why high 
altitudes are reached by this rocket dust storm. Con- 
versely, we can simulate with the LMD-MMM what 
would have happened if the OMEGA storm had ap- 
peared in dusty northern winter, e.g. at L s = 240°. 
This solar longitude is chosen so that background 
dustiness is significantly increased compared to the 
reference case (rb ac k = 0.6 — 0.7 instead of 0.3 — 0.4), 
while incoming solar energy is approximately sim- 
ilar. Results are shown in Figure [13] The dust 
disturbance does develop a convective behavior, but 
the warmer environment causes smaller values of 
CAPE C and vertical acceleration than in the refer- 
ence case. Hence the rocket dust storm reaches lower 
altitudes (15 — 20 km) and loses its vertical extent 
earlier in the afternoon. The clearest seasons with 
respect to background dustiness appear more propi- 
tious to high-altitude rocket dust storms. 

Two additional factors, while not as crucial as 
the two previous ones, could influence the develop- 
ment of rocket dust storms: atmospheric lapse rate 
and wind shear. Those are related to vertical vari- 
ations of (respectively) temperature and horizontal 
winds. A more stable profile would tend to in- 
hibit convective ascent of rocket dust storms. As 
the Martian atmosphere lacks a stratosphere, rocket 
dust storms are probably devoid of anvils at their 
summits that are common in terrestrial cumulonim- 
bus. Yet contrasts in atmospheric stability could 
cause spatial and temporal variability in the alti- 
tude reached by rocket dust storms. This proba- 
bly combines to the effect of background dustiness 
to make northern winter the least propitious season 
to high-altitude rocket dust storms (unless exception- 
ally dusty). Dustier background atmosphere implies 
warming in mid-troposphere and cooling in lower tro- 
posphere, hence a more stable temperature profile 
which would tend to mitigate vertical convective ac- 
celeration (as is the case in Figure 13). As far as 
horizontal wind shear is concerned, strong vertical 
shear could make rocket dust storms lose their ver- 



tical extent, but this does not seem to happen when 
convective acceleration is significant. For instance, a 
strong near-surface vertical shear of horizontal wind 
does not prevent the OMEGA storm to develop and 
rise (see Figures [4] and 10). 

We conclude that rocket dust storms would reach 
particularly high altitudes in tropical regions in 
northern spring and summer, where and when at- 
mospheric dustiness is low. Detached layers of dust 
which result from the ascent and horizontal spread- 
ing of rocket dust storms, as is the case for the 
OMEGA storm (section EJ), would tend to follow 
the same spatial and seasonal variability. In other 
words, if we assume that rocket dust storms are the 
source of detached layers of dust, we expect the lat- 
ter to be particularly discernible in tropical regions 
in norther spring and summer. This is in agree- 



ment with the conclusions of Heavens et al. (2011a 



based on MCS measurements about the presence of 
a so-called "high-altitude tropical dust maximum" 
throughout northern spring and summer. This max- 
imum has a greater magnitude and altitude in the 
northern hemisphere. This can be easily understood 
since in spring/summer the northern hemisphere re- 
ceives more incoming sunlight than the southern 
hemisphere, hence is more prone to strong rocket dust 
storms conducive to high-altitude detached layers. 



Furthermore, Heavens et al. (2011a) found in MCS 



measurements that the tropical dust maximum is par- 
ticularly well-defined in the range L s = 110 — 160° 
for MY28 and L s = 45 - 140° for MY29, and related 
this difference to THEMIS and MARCI observations 
which showed a lack of "early season" tropical dust 
storm activity in MY28 compared to MY29. Finally, 



Heavens et al. (2011b) reported the significant longi- 



tudinal variability of the high-altitude tropical dust 
maximum, which maintained over time scales so large 
that dust lifting, transport and removal processes (i.e. 
fast mesoscale phenomena) must oppose advection, 
sedimentation and eddy diffusion (i.e. slow large- 
scale phenomena) . This is in line with detached layers 
of dust being related to rocket dust storms. 

Limb sounding by TES (iClancy et all 120101 e.g. 



their Figure 10) and MCS ( |Heavens et al.| |2011a| 



e.g. their Figure 5) show that northern winter (dusty 
season) is not devoid of detached layers of dust. The 
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Figure 13: Same as Figure [4] except LMD-MMM simulation is performed at L s = 240°. 
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main distinction with northern spring/summer is how 
the altitude of enriched layers of dust compares with 
the summit of the background dust layer (roughly 
equivalent to what Heavens et al. (2011a) named 



falloff height). In northern spring/summer in tropi- 
cal regions, the summit of the background dust layer 
is 40 km while it is 60 — 80 km in winter (e.g. Fig- 
ure 4 and 5 in Montmessin et al. 2006). Hence en- 



riched layers of dust are detected within the back- 
ground dust layer in northern winter while they ap- 
pear close to the summit of the background dust layer 
in spring/summer (if not above it). This seasonal 
contrast in the nature of detached layers is best sum- 



marized by Figure 1 in Heavens et al. (2011b). This 



difference can be related to the fact that background 
dustiness is larger in northern winter, which causes 
rocket dust storms to reach lower altitudes compared 
to northern spring/summer (Figure 13). In other 



words, rocket dust storms would tend to create en- 
riched layers of dust within background dust layer in 
dusty season; and to form "truly detached" layers of 
dust in clear season, such as the high-altitude tropical 
dust maximum witnessed by MCS. 



nitude to first order (section 4.1), but could also ac- 



count for their observed diurnal and seasonal variabil- 
ity, and to some extent their interannual and longitu- 
dinal contrasts (section |5.1| . This posits that rocket 
dust storms are plausibly the best explanation for 
detached layers of dust in the Martian atmosphere. 

One of the difficulties found by |Heavens et al.| 
(2011b) with this scenario is that MOC statistics of 



dust storm activity (Cantor et al. 2001 2006) show 



that tropical dust storm activity is maybe too low 
to support the high-altitude tropical dust maximum. 
This difficulty is not irreconcilable with our conclu- 
sion that dust storms are a major source of detached 
layers of dust. 

1. The limited horizontal extent and extremely 
fast evolution of rocket dust storms could pre- 
vent their detection by remote-sensing nadir 



techniques with large footprints (THEMIS, TES) 
which can only monitor background dustiness 
and not mesoscale contrasts. We also found 
rocket dust storms possibly start their ascent 
in late morning instead of early afternoon (sec- 
and 



tions 



4.2.2 



4.2.3), which could prevent 



MGS/MOC, operating at local time 1400, from 
detecting most of those storms. 

2. Not a large number of rocket dust storms are 
required for detached layers to be significantly 
loaded with dust particles. Our simulations show 
that a single rocket dust storm (constrained 
through observations) yields within a few hours 
a detached layer of dust which corresponds to 
typical values of density-scaled optical depth ob- 
served by MCS. Furthermore, we showed how 
detached layers of dust could survive for several 
sols once formed by rocket dust storms owing to 
slow sedimentation and daytime radiative warm- 
ing of dust particles. Even intermittent rocket 
dust storms could lead to high-altitude detached 
layers of dust observable by MCS throughout 
northern spring/summer. 



5.2 Sources of detached layers of dust 3 - 

Not only rocket dust storms have the ability to form 
detached layers of dust at correct altitude and mag- 



MOC statistics might be incomplete: tropical 
dust storms were observed in the tropics by 
Viking (cf. Figure 13 by |Cantor et al.[ |2001[); 



the OMEGA storm itself was detected in an area 
where dust activity is particularly low according 
to MOC, as is also the case for the convective- 
looking dust storm caught by the HiRISE cam- 



era in Figure 9 of Malin et al. ( 2008 ); dust optical 



depth measurements on board Mars Exploration 
Rovers in low latitudes exhibit sudden large vari- 
ations in northern summer which are reminiscent 



of dust storm activity (Figure 13 in Smith et al. 
20061). 



Other potential sources for MCS detached layers 



of dust were explored by Heavens et al. (2011b) 



Mesoscale modeling by Rafkin et al. (2002) showed 
that anabatic (upslope) winds converging at the top 
of Arsia Mons lead to transport of dust particles at al- 
titudes 20 km over large horizontal distances. While 



their contribution cannot be ruled out, Heavens et al. 



(2011b) noted that those cannot account for detached 
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layers of dust in areas distant from significant topog- 
raphy. A distinct scenario is that mesoscale circula- 
tions are not even involved in the formation of de- 
tached layers of dust: scavenging by water ice clouds 
could be a plausible explanation. Large water ice col- 
umn optical depths, forming the aphelion cloud belt, 
are indeed observed when the high-altitude tropical 



dust maximum occurs (Wang and Ingersoll 2002 



Clancy et al. 2003). Heavens et al. (2011a) however 
resolved tropical maxima with MCS when the aphe- 
lion cloud belt is dissipated. We complement this 
statement by noting enriched layers within the dust 



background layer in northern winter (Clancy et al. 



2010 ) cannot be accounted for by scavenging through 
water-ice clouds which are inexistent in the lower 
troposphere at this season. Heavens et al. (2011b) 



and |Rafkin| (|2012]) also argued that scavenging alone 
could not yield the observed high-altitude detached 
layers of dust with mixing ratios larger than in the 
lower troposphere. Significant mesoscale transport 
of dust particles is necessary for this to occur. Our 
mesoscale simulations of rocket dust storms are in 
line with those conclusions. Scavenging by water-ice 
clouds would only help to reinforce an existing de- 
tached layer of dust. 

To explain the occurrence of high-altitude tropi- 



cal dust maximum, Heavens et al. (2011b) also pro- 



posed to extend the idea of Fuerstenau ( 2006 1 that 
solar warming of suspended particles in dust devils 
could lead to convective acceleration and dust parti- 
cles reaching altitudes as high as 15 — 20 km. This 
scenario suffers several possible limitations. 



1. 



Compared to our computations for a t ~ 4 
storm, heating rates estimated by Fuerstenau 
( |2006[ ) within dust devils appear overestimated, 
especially given that dust devils are of lower op- 
tical depth than local dust storms and that large 
dust loading is mostly found at the edges of the 
vortex. In addition to this, THEMIS images in 



Towner ( 2009 ) indicate that dust devil cores are 



mostly colder than their environment. 

2. While dust devils imply upward transport of 
large quantities of dust particles, an unrealisti- 
cally large number of them would be necessary 
to advect as many dust particles as a rocket dust 



3. 



storm extended over tens of kilometers such as 
e.g. the OMEGA storm (Figure [IJ. 

Should dust devils be involved, regions propi- 
tious to their apparition (e.g. Amazonis Plani- 
tia, Fisher et al. 2005) would play a role in 



the spatial variability of detached layers of dust, 
which is not clearly observed. 

4. Considering dust devils as individual convective 



towers as done in Fuerstenau ( 2006 ) and Heavens 



et al. (2011b I is not an entirely correct approx- 
imation. Dust devils form when dust get trans- 
ported within convective vortices which arise as 
part of the complex turbulent growth of the un- 
stable daytime PBL (|Kanak et aT||2000|). The 



adaptation of the simplified convection model of 



Gregory (2001) to Mars proposed by Heavens 
et al. (2011b) would actually be more applicable 
to rocket dust storms than to dust devils. This 
is acknowledged by the authors themselves, as 
well as the need to address questions related to 
dust devils through turbulent-resolving simula- 
tions that include radiatively-active transported 
dust particles. 



5.3 Implications of dusty deep convec- 
tion on Mars 

Deep convection triggered by radiative effect of dust 
particles is a remarkable feature of the Martian me- 
teorology with numerous implications. 

Rocket dust storms, and their plausible intimate 
link with detached layers of dust, underline the im- 
portance of mesoscale phenomena in driving the dust 
cycle on Mars. This does not contradict but comple- 
ments the role played by planetary-scale circulations 



in dust transport ( 


Newman et al. 


2002 Kahre et al. 


2006 


Madeleine et al. 


2011 


1. Our mesoscale simu- 



lations show indeed how detached layers of dust re- 
sult from both rocket dust storms (for vertical trans- 
port) and large-scale horizontal winds (for horizon- 
tal spreading of dust particles). A notable conse- 
quence is that GCMs lack both temporal and spa- 
tial resolutions to simulate local storms, hence can- 
not accurately represent the vertical transport of dust 
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particles in the troposphere. Heavens et al. (2011a 



noted that the observed dust distribution at north- 
ern summer, notably the high-altitude tropical dust 
maximum, cannot be reproduced by GCM simula- 



tions of Richardson and Wilson (2002) and Kahre 



et al. (2006), while agreement is better in northern 



winter. This can be understood as rocket dust storms 
are likely to be more active and reaching higher al- 
titudes in the former than in the latter season (see- 
In a more general sense, quite similarly 



5.1) 



tion 

to moist convective cells in terrestrial climate, rocket 
dust storms could significantly contribute to vertical 
transport in the ascending branch of the Hadley cell 



(Rafkin 2012) 



The analogy with moist convective storms on Earth 
opens broader perspectives. We describe in section [4] 
how radiatively-induced deep convection in the Mar- 
tian atmosphere significantly impacts wind and tem- 
perature fields. Rocket dust storms, as their terres- 
trial counterparts, are efficient ways to transport heat 
and momentum. We also show that detached layers 
of dust could be a significant source for atmospheric 
warming and supply of convective energy several sols 
after being formed by a rocket dust storm. Thus we 
expect rocket dust storms to have a strong impact on 
the Martian climate which will need to be taken into 
account in GCMs. This applies to 



thermal structure: Madeleine et al. (2011) no- 



ticed their GCM underestimates temperature at 
altitudes 20 — 40 km during rapid increases of 
column dust optical depth; 



A parameterization of mesoscale dust storms, simi- 
lar to cumulus parameterization on Earth, appears 
necessary to improve GCM predictions on Mars. 

The local impact of rocket dust storms is as impor- 
tant as their influence on the global climate. A po- 
tentially crucial element is the emission of mesoscale 
gravity waves by rocket dust storms. Gravity waves 



induce large temperature perturbations (Spiga et al 



2012) and, especially when they break, exert a drag 



on large-scale circulation which needs to be pa- 
rameterized in Martian GCMs (Forget et al. 1999[ 



Medvedev et al. 2011 ). By accounting for the convec- 



tive source in addition to the (commonly considered) 
topographical one (Creasey et al. 2006), mesoscale 



simulations of rocket dust storms could help to im- 
prove both the interpretation of GW-induced phe- 
nomena, and GCM parameterizations. Mesoscale dy- 
namical disturbances implied by rocket dust storms 
are actually crucial for the exploration of Mars in gen- 
eral. Figures[5]and[7]showed that large horizontal and 
vertical winds are associated with the occurrence of 
rocket dust storms. Modeled winds reach limits over 
which entry, descent and landing might not be safe 
(Vasavada et al. 2012 1. Thus the possibility that 



deep convection occurs in the Martian atmosphere 
shall be considered for designing and engineering fu- 
ture missions in the Martian environment. Rocket 
dust storms are all the more critical for the explo- 
ration of Mars as those could produce strong electric 
fields through triboelectric charging and give rise to 



lightning and radio emission (IRenno et al. 2003 Ruf 



et al. 2009 ), as is the case for thunderstorms on Earth 



and other planetary atmospheres (Russell et al. 2011 



predictability: Rogberg et al. (2010) showed how Fischer et al. 2011 ). 



meteorological predictability error could grow 
large as a result of uncertainties in the dust dis- 
tribution; 

• atmospheric dynamics: vertical transport of heat 
and momentum in rocket dust storms could per- 
turb thermal tides ( Lewis and Read 2003[ ) , trig- 



ger Rossby waves (especially within tropical lat- 
itudes, e.g. section 2 in Schneider and Liu 2009 ) 



and couple to large-scale circulation to give rise 
to instabilities (e.g. the Madden-Julian Oscilla- 
tion on Earth, Zhang 2005). 



Rocket dust storms do not only impact dust cy- 
cle and atmospheric circulation. Strong convective 
motions would induce fast upward transport of any 
molecule or aerosol present in the vicinity of the dust 
storm (before high-altitude jet streams induce exten- 
sive horizontal transport) . A notable consequence for 
the Martian water cycle is that rocket dust storms 
could cause deep and efficient vertical transport of 
water vapor. This might help to explain the (of- 
ten rapid) surge of detached layers of water vapor 
(Maltagliati et al., submitted to Icarus) and water ice 
clouds ( Clancy et al. , 2009 Vincendon et al. 2011 1 at 



2G 



altitudes above 40 km. Possible upward transport of 
molecular species and aerosols by rocket dust storms 
has numerous implications too for photochemical cy- 



cles on Mars (Rafkin 12012), especially if electrical 
activity is significant (|Farrell et al. 20061. Rocket 



dust storms would also control the amount of conden- 
sation nuclei in the atmosphere, which is key to un- 



derstand supersaturation on Mars ( Maltagliati et al 



2011) and heterogeneous chemistry (Lefevre et al. 
2008). 



The convective behavior evidenced here for local 
dust storms on Mars also apply to larger dust storms. 
Dust fronts, which develop at the edge of polar caps 
and as a result of baroclinic instability, sometimes 



exhibit convective-like appearance ( Strausberg et al 



2005). Regional dust storms might be for rocket dust 



storms what mesoscale convective systems are for in- 



dividual thunderstorms on Earth (Houze 2004 1. A 



behavior similar to rocket dust storms, or a cluster 
thereof, might be involved in the development of Mar- 
tian global dust storms, whose mechanisms are left 
to be explained. In addition to other factors such 



as the availability of dust reservoirs (Pankine and 



Ingersoll 2002), the role played by mesoscale phe- 



nomena might explain why GCMs have difficulties re- 
producing the intcrannual variability of those events 
(iNewman et all 120021 iBasu et all 120061). This is 



a scientific problem probably worth being explored 
in the future with mesoscale models which, contrary 
to GCMs, have the ability to resolve convective mo- 
tions. 
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